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Abstract 

The dynamics of decisions in complex networks is studied within a Markov process frame- 
work using numerical simulations combined with mathematical insight into the process mech- 
i-C ' anisms. A mathematical discrete-time model is derived based on a set of basic assumptions 

on the convincing mechanisms associated to two opinions. The model is analyzed with respect 
to multiplicity of critical points, illustrating in this way the main behavior to be expected in 
the network. Particular interest is focussed on the effect of social network and exogenous mass 
media-based influences on the decision behavior. A set of numerical simulation results is pro- 
vided illustrating how these mechanisms impact the final decision results. The analysis reveals 
(i) the presence of fixed-point multiplicity (with a maximum of four different fixed points). 
C/2 ■ multistability, and sensitivity with respect to process parameters, and (ii) that mass media 

^~>| have a strong impact on the decision behavior. 



Or 



1 Introduction 

Recently, the problem of modeling, analysis and simulation of rumor and decision dynamics in 
complex networks has obtained increasing attention, including studies on rumor spread in social 
and small- world networks [T], and decision dynamics in scale- free networks [2]. While the problem 
of rumor spreading [T] considers the time evolution of rumor spreaders and stiflers, the decision 
dynamics [2] is concerned with the time evolution of different decisions in terms of the associated 
competitive interplay, similar to the one observed in biological studies on different predators com- 
peting for the same prey, or chemical species competing for the same reactant. The most notable 
difference with respect to these classic fields of competition dynamics consists in the fact that in 
social system dynamics the topology of the underlying contact network is substantial. For example, 
social phenomena involve different personalities with different number of contacts each, so that 
the dynamics take place with an intrinsic distributed character and, in particular, emergence of 
nearly homogeneous groups (so called clusters) is possible. These differences make it necessary to 
analyze social behavior from the view point of dynamic networks. It is particularly noteworthy 
that it has been observed that the underlying network topology, unless changing in time, always 
fulfills the same characteristics in terms of node degree distribution (see e.g. [31 2J [5]). For the 
purpose of understanding the mechanisms underlying dynamic phenomena over dynamic networks, 
mathematical models can be developed and analyzed formally as well as with numerical simulation 
studies. 
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With particular emphasis on the dynamics of decision competition in scale-free networks, in [5] 
the regularity of spreading of information and public opinions towards two competing products in 
complex networks was analyzed. They proposed a simple linear model and simulated the associated 
decision trajectories over time. It was highlighted that, in contrast to most models studied so far via 
modified SIS and SIR models [HI El EH HH US HU , in a real life setting there are frequently 
various information being spread simultaneously (diverse virus, multiple opinions, rumors, etc.), 
and, in contrast to the single opinion spread dynamics, these may be mutually strenghtened or 
even annihilated. The competition-based dynamics of two opinions which freely flow in a complex 
network was studied, showing some interesting and important facts concerning the behavior of de- 
cision competition. Nevertheless, studying the underlying mechanisms of mutual influence between 
nodes in the network is expected to lead to a nonlinear dynamics, just as in the SIS and SIR case 
mentioned above. Accordingly, it may be supposed that the inherent competition mechanisms lead 
to classical nonlinear phenomena such as multiple attractors and parameter sensitivity. 

In the present study a mathematical model is introduced and analyzed for the dynamics of 
two competing opinions (with neutral intermediate state) which explicitly accounts for the kind of 
nonlinear interactions inherent to the dynamics of competing opinions in complex social networks. 
A generalized model for the interaction mechanisms is employed which has been recently proposed 
in [13j and provides intermediate cases between the classical contact and reactive process, and the 
influence of exogenous mass media opinion propagation into the network is considered explicitly. 
The mathematical model is analyzed with respect to attractor multiplicity, delimiting in this way 
some system inherent behavior possibilities. In particular the study focusses on the effect of opinion 
adaptation through (i) social contacts, and (ii) exogenous mass media influence. In order to analyze 
these effects, a series of numerical simulations is provided for the- two extreme case's of contact and 
reactive process [13] . illustrating typical behavior in presence of: 

• static contacts 

• dynamic contacts 

• exogenous mass media propagation of one opinion 

• exogenous mass media propagation of both opinions. 

The presented results illustrate (i) typical nonlinear behavior such as attractor multiplicity, and 
sensitivity with respect to process parameters, and (ii) the strong impact mass media have on the 
decision behavior. The presented results establish an extension of those reported in [2j. 

From a methodological point of view, we are putting together the quantitative and qualitative 
potential of mathematical modeling and simulation, with some basic concepts of dynamical systems 
theory on attractor multiplicity and stability on the basis of numerical simulation studies. 

The paper is organized as follows. In Section 2, the nonlinear opinion dynamics model for 
dynamic networks is derived based on a set of basic assumptions, and the fixed points are established 
for the corresponding static and unperturbed version. In Section 3, numerical simulations are 
presented which illustrate the main behavior observed in decision dynamics in social networks. In 
Section 4, the main contributions of the study are summarized. 
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2 The decision dynamics model 



In this section a mathematical model for the prevalence of two opinions, A and B, with a neutral 
intermediate state N is derived based on the Markov process assumption [T5] that the state at any 
time instant t > 1 (t £ N) depends exclusively on the state at the immediately preceding time 
instant t— 1. The result is a discrete time 3-/V-dimensional nonlinear model, where N is the number 
of nodes in the network. 

2.1 Basic assumptions and nomenclature 

In order to derive the model we introduce the following basic assumptions: 

(Al) Two opinions A and B arc being propagated in a social network with N nodes and Power 
Law degree (k) distribution [TB] 

P[k] ~ fc~ 7 (1) 

(see Figure [1]). The corresponding network topology is reflected in the square adjacency 
matrix A(t) with time- varying entries Aij(t) £ {0, 1}. 




-2.63 



Figure 1: Power law degree distributed network with N = 100 nodes, and P[k] < 
(A2) The total number of nodes N is constant. 

(A3) Any node i can be in one of three states: A if it has the opinion A, B if it has the opinion B, 
and Af if it is neutral with respect to the opinions A and B. 
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(A4) In order to change from opinion A to B (or vice versa), the node has to pass through the 
neutral state AT, i.e. there is no direct connection between the states A and B. The associated 
state-diagram is depicted in Figure [2] 

(A5) The probability for a node i to be in the state A, B, or TV, is denoted by flj, bi and n,, 
respectively. It holds that aj, bi, rii € [0, 1] for any i = 1, . . . , N. 

( A6) If the node i has the opinion A (or B) , it will convince of opinion A (or B) the nodes connected 
with him with probability ka (or Kb)- 

(A7) A neutral node i (i.e. a node which is in the state A/"), does not convince any neighbor node. 

(A8) The impact any neighbor node j ^= i has on node i is independent of j. 

(A9) The influences of any two nodes j ^ k =/= i on node i are mutually independent. The total 
effect of all neighbors of node i is given by the mean influence. 

(A10) The opinion A (or B) is interchanged between two nodes i and j Xa (or Xb) times in each 
time step [15] . 

(All) The mass media influence can be modeled in terms of an exogenous perturbation with period 
of appearance T, taking into account the underlying automata transition possibilities (see 
assumption A4 and Figure [21). 




Figure 2: State transition diagram. 
A direct consequence of assumptions A2 and A4 is that 

cti(t) + h(t) + m(t) = l. (2) 

Assumption A8 is motivated by the fact that any node j has the same impact on i (assumption A7) 
and the final state of node i after contacting his neighbors will be a weighted sum of all particular 
contacts. Assumption A9 is associated to the type of interchange mechanism |13j . Two extreme 
cases have been reported in the literature: the contact process (with Afc = 1, k = A, B), and the 
reactive process (with X k —> oo, k = A, B). In [T3] a generalized model for the contact rate was 
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proposed covering all intermediate scenarios. In the present study the model proposed in |13j for 
the contact rate is adapted by defining the transmission rates associated to opinion A and B as 

where iVj(i) is the total number of neighbours of node i at time i, i,e. 

N i (t) = Y^A itj (t). (4) 

For the reactive process (Afc — > oo, k = A,B) = 1, k = A, B. The transition probabilities 
associated to the process are denoted as follows: 

• If the node i is in the state A, B or Af, the probability for node i of remaining in A, B, or Af 
is denoted by on, fti, or Vi, respectively. 

• If the node i is in the neutral state Af, the probability of being convinced of opinion A (or B) 
is denoted by fa (or <x;). 

2.2 Derivation of the Markov process model 

Having as point of departure the preceding assumptions, the transition probabilities on, fa 
and <Ji are determined. 

As representative case, suppose that a node i is of opinion A and has contact with a node j =/= i. 
The following events are possible: 

• The node j is in the state A (or AT) with probability a,j (or tij), and consequently will not 
convince node i of changing its opinion. 

• The node j is in the state B with probability bj, and consequently is trying to convince node 
i with probability of success given by t^kb, and will fail to convince i of opinion B with 
probability 1 — rf^KB, with rf^ given in ©. 

On the basis of these considerations, the probability for node i of not being convinced by node j of 
changing from state A to Af (the only alternative to remaining in A) , is given by the sum 

ctij = cij + n,j + (1 - r?K B )bj = 1 - r^K B bj, (5) 

where the last equality is a consequence of property (|2|). 

The same reasoning applies to the probability /3y of node i being in state A of not being 
convinced to change from the state B to the state Af by its neighbor j / i [having in mind property 
©]: 

N = b j + n j + C 1 - rfjKAjaj = 1 - r^K A a v (6) 

In the case that node i is in the state Af and in contact with a node j ^ i, the following events 
imply that i remains in Af: 
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• The node j is of opinion M with probability rij, and does not convince i of changing its 
opinion. 

• The node j is of opinion A (or £>) with probability a,j (or bj, and does convince i of opinion 
A (or B) with probability r^KA (or rfjUs), and fail with probability \ — t^ka (or 1 — rf^Ks)- 

Accordingly, having in mind property ([2]), the probability for node i of not being convinced by node 
j ^ i is given by 

Vij = rij + (1 - r^KA)a,j + (1 - rfjKsJbj = 1 - r^KA^ - r^Ksbj, (7) 
and the probabilities of node j ^ i convincing node i of opinion A or opinion i? are given respectively 

by 

Hij =rfjKAaj, <Xij =rfjK B bj. (8) 

Next, to take into account the combined effect of all neighbors of node i on the transition prob- 
abilities, recall assumption A8. The overall transition probabilities for node i at time t are given 
by: 

A(*) = j^E^Mt) = ^yE^Ji - ^(t^WI, 

= ^ E 3 -* f«W = E^J 1 - r^(i)KAaj(t) - rf 3 (t) KB b 3 (t)}, (9) 

= -^E,W*) = ^E^W^M- 

The associated discrete time model for the dynamics of decisions between .4 and 23 in the network 
is given by 



Oj(i + 1) = ai(t)ai(t) + Hi(t)rii(t), a;(0) = a i0 

TH(* + 1) = i^H(t) + [1 - <*(<)]<*(*) + [1 - ft(f)M*). «i(0) = n, l0 

h(t + l) = pi{t)k{t) + <n(t)m(t), bi(0) = b i0 
o = i- [ai{t) + bi(t) + m(t)]. 



(10) 



Taking into account the identities [sec Q] 

fii(t) = 1 - &(*), a i (t) = l-a i (i) (11) 

model ([TU)) is equivalent to 

t»i(t + 1) = Oi(t)ai(t) + [1 - ft(t)][l - a»(*) - h(t)], ai(0) = a, 
6 i (t+l)= i 8 i (t)l> i (t) + [l-a i (t)][l-a i (t)-6 i (t)], 6,(0) = 6 i0 (12) 
n ( (t) = l-o ( (t)-6 1 (t) 

with aj(t) and /?i(t) defined in (|9]). In vector notation the preceding dynamics are written as 

Xi(t + 1) = f[xi(t)], Xi(0) = x M , Xi = [a i ,b i ,n i ]€T, i = l,...,N, (13) 
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where T is the two-dimensional triangle set (see Figure 

T = {z= [z l7 z 2 ,z 3 }' £ [0, l] 3 C R 3 | zi + z 2 ,+z 3 = 1}. (14) 
Introduce the mean decision probabilities 

^ A/" ^ A/" j JV 

i—1 i—1 i—1 

and the associated mean probability p(t) = [p a (t), Pb(t), p n (t)]' , it can be easily seen that 

p{t)eT, Vt>0. (16) 

In order to analyze the influence of mass media on the decision behavior, consider the time- 
varying index subset 

4(t)c{l,...,4 #Ik = [VkX\, k = A,B, (17) 

affecting rj% of the total population (here indicates the cardinality of the set /) , and introduce 
the associated characteristic function 

**<*>={!; Wm ' k = AB - (18) 

In terms of the index subset Ik, k = A,B the mass media influence can be modeled as follows: 



ai {t + 1) = ai(t)oi(t) + [1 - ft(t)][l - o,-(i) - 6,-(t)]+ 

(tK(i)n*(i) - Xf (t)v B (t)oi(t), Oi(p) = a M 
bi(t + 1) = f3 l (t)b l {t) + [1 - Oi(t)][l - cn(t) - 6i(*)]+ 

+Xf (*)« fl (t)ni(t) - xt(t)v A (t)h(t), 6,(0) - 6 i0 



(19) 



where 



is a function which is nonzero only in discrete time instants sXfc, s £ N, k = A, B with period T^. 

It should be noted that, in comparison to the study presented in [2], there are some substantial 
differences to the present work. First, the model built in [3J is linear while the model (|T5)) is 
intrinsically nonlinear. Second, the model (|T9l allows to take into account different modalities of 
transmission process between the two extreme cases of contact process (A = 1) and reactive process 
(A — ► oo ) as proposed in [13]. Third, in (|1Q|) the influence of exogenous opinion propagation into 
the network is explicitly considered. 



2.3 Basic characterization of the dynamic behavior 

Assume for the moment that the connections between individuals do not vary over time (i.e., the 
adjacency matrix A is constant, and hence Ni is constant for any node i). Based on this assumption 
the fixed points associated to the dynamics ([TO]) can de determined in order to establish limit case 
conditions for the time evolution of the opinions A and B in the network. 
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The fixed point condition read 

a.i(t + 1) = Oi(i) = at, bi(t + 1) = bi(t) = bi, n^t + 1) = rii(t) = Hi 



(21) 



where <Xi, 6j, are the fixed point values of the stochastic variables a,i(t), bi(t) and rii, respectively. 
It follows that 



= (oti - l)a t + (1 - fiijUi, 
= (A - 1)6* + (1 - a*H, 
= -(1 - a*H - (1 - Pi)lH + (2 - a t - fam 



(22) 



Given that for a, = (or &j = 0) for all i it holds that fa = 1 (or a, = 1), it follows that the three 
vertexes 



(23) 



of the triangle set T (| 14[) (Figure [3]) are particular fixed points. Furthermore, solving the equation 
set (|22|) in general with respect to di, bi and m yields the following set of implicit solutions 





' 1 " 




' " 




' " 


X = 





, x = 


1 


, x = 

















1 



1 



-, bi 



"""tf + fc + i' 1 g + 6 + 

parameterized by the scalar 



rii = 1 - en - h 



tf + fc + l 



, i = l,...,JV, (24) 



1 j¥=i 1 i¥* 



Introducing the solution in vector notation 

Si = [m, bi, rii}' 6 T, 
the following limit cases can be directly derived: 



ka — > 0, => £j — > oo and Sj 



, kb — > => £, — > and 



(25) 



(26) 



(27) 



corresponding to the two bottom vertexes of the triangle set T (fT4|) (Figure [3]). In terms of the 
opinion density p (|16|) there are four associated fixed points 
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Ei=i ^ 



(28) 



where ai,bi and n, are given in (|24[) . Naturally, the limit cases (|27|) can be expressed in terms of 
the mean probability solution s p as 



KA ->■ 0, 







TV ^ 



(29) 
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Figure 3: Sketch of the triangle set T (fT4)) where the decision probabilities of each node i and the 
mean probabilities p evolve, and the curve S parameterized by the scalar £ ([251) . 

Summarizing, for any parameter combination (ka, Kb, Xa, Xb) > there are four fixed points, 
and in the limit ka = (or kb = 0) there are three fixed points. This fact illustrates the multiplicity 
of fixed points and sensitivity with respect to the dynamic's parameters, and reveals that the 
nonlinearity introduced by the interchange mechanisms is too strong to be neglected. This fact 
establishes a main difference with respect to the recent study [2] , where a linear model was used for 
approximating the decision dynamics. A formal analysis of the stability properties of the four fixed 
points nevertheless goes beyond the scope of the present study, and here we circumscribe ourselves 
to numerical simulation studies in order to find the main behavior expected in decision dynamics 
in social networks, and study the impact of mass media exogenous perturbations. 

3 Simulation results 

In this section, the dynamic behavior of decisions in social networks is analyzed for four different 
scenarios: 

• Static network, i.e. static adjacency matrix with power law (scale free) distribution. 

• Dynamic network, i.e. time-varying adjacency matrix with power law distribution over any 
time interval. 

• Dynamic network with exogenous mass media influence of one opinion on 18 % of the popu- 
lation with period Ta = 8 time units. 

• Dynamic network with exogenous mass media influence of both opinion on 10 % of the pop- 
ulation with period Ta = 8 time units and Tb = 16 time units. 
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The simulations were taken out for a total population of N = 10.000 nodes, and a power law node 
degree k distribution ([T]) with [16] 



7 = 2.16. (30) 
3.1 Static network without exogenous perturbations 

In order to illustrate the main (nonlinear) behavior of the decision dynamics process according to the 
fixed point multiplicity discussed in Section 2.3, the two limit cases with A^. = 1, k = A, B (contact 
process) and — s- oo, k = A,B (reactive process) arc illustrated for three different parameter 
scenarios. For the purpose at hand the projection of the trajectories in the triangle set T (|14[) 
onto the (p a , p;,)-plane is presented. Note that accordingly it is possible that trajectories seem to 
intersect in the projection. 

For the purpose at hand consider the following three parameter sets: 

ka = 0.1, Kb = 0.9 (ka < kb), ka — kb — 0.5 (ka = Kb), ka — 0.9, kb — 0.1 (ka > Kb) (31) 

and that A and B are propagated via contact processes (i.e., Xa = Ab = 1)- The corresponding 
simulation results are presented in Figure [H 




Figure 4: Projection onto the (p a , pb)-plane of the trajectories of the mean decision probability p 
(|16p for the three different parameter sets given in pip and a contact transmission process. 

It can be seen in Figure U] that: 

(i) for the first parameter set (ka < Kb) ah the trajectories move towards an attractor in the 
upper left corner of the (pA, ps)-plane. 

(ii) for the second parameter set (ka = Kb) there are two attractors (one in the upper left and one 
in the lower right corner) towards which the decision trajectory p may converge, in dependence 
on the initial value. 

(iii) for the third parameter set (ka > Kb) ah the trajectories move towards an attractor in the 
lower right corner of the (pa, ps)-plane. 
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According to these results, there is a strong sensitivity, in the kind of structural instability |17j . of 
the decision behavior in the network on the convincing parameter pair (ka, Kb), manifesting itself 
in particular through the fact that, in the passage from the first to third parameter set in pip, the 
attractor at the lower right corner initially is unique, than coexists with a second attractor in the 
upper left corner, and finally dissapears leaving the attractor in the upper left corner as the unique 
one. Hence, for some parameter combinations the final decision state will depend strongly on the 
initial condition. 

It should be noticed that close to the fixed-point, the dynamics becomes very slow in corre- 
spondence to the fact that almost the whole network is of one single opinion, so that the existing 
contacts only confirm the present decision state. This behavior is predicted by the model (fT2"j) , 
given that close to the fixed point the functions aj,/3j ~ 1 for any node i, and hence the system 
attains a slowly varying (i.e., quasi steady-state) behavior according to 

Oj(t + l) (*ai(t), bi(t + l)^bi(t), ni(t + l)«fii(t). (32) 

Our conjecture is that this fact is also related to the generation of decision clusters within the 
network. 

In Figure [S] are presented the simulation results for the three parameter combinations defined 
in (j3"T|) . but for reactive processes with = Xb — > oo, or equivalently rf- = rfj = 1. It can 
be observed that the global over-all behavior is quite similar to the one observed for the contact 
process. Comparing the behavior over time, the only substantial difference consists in that the 
convergence speed is considerably faster for the reactive process, in accordance to the fact that the 
transmission rate is much higher. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 

P. P. P. 



K A < K B K A — K B K A > Kg 

Figure 5: Proyection onto the (p a , ph)-plane of the trajectories of the mean decision probability p 
(fro]) for the three different parameter sets given in (f3"Tj) and a reactive transmission process. 

3.2 Dynamic network without exogenous perturbations 

Next, consider the case that the network connections vary with time in such a way that the over- 
all degree distribution ([1]) is maintained. In this case the system dynamics (|12|) are no longer 
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autonomous, and thus no fixed points exist for the system (unless in the mean, i.e. in terms of 
Pk,k = A, B, N, asymptotical stability-like behavior may be observed). 

The associated decision time evolution is compared in Figure |6] to the static network case for 
the parameter value and initial condition set 

ka = 0.1, =0.9, Pa (0) = p b (0) =0.5, |O„(0) = (33) 

and a reactive process (i.e., = l,k = A,B). It can be observed that the over-all behavior is 
quite similar, with the important difference that for the dynamic network case the trajectories are 
smoother and converge in about 9 time steps, while in the static network set-up still the fixed-point 
is not reached. These observations are in line with the results presented in |18j . and are probably 
a result of the fact that in the dynamic network set-up, i.e. with contacts changing at each time 
instant, the formation of decision clusters is less probable than in the static network set-up. This 
conjecture should be carefully studied in future research. 




5 10 15 20 5 10 15 20 

Time t Time t 



Static Network Dynamic Network 



Figure 6: Time evolution of the mean decision trajectories associated to decision A (thin grey 
line) and B (thick black line) for the parameter value and initial condition set (|33j) and a reactive 
transmission process in a scale-free dynamic network. 

3.3 Dynamic network with mass media influence of one single opinion 

To analyze the impact of mass media on the decision behavior of a dynamic network, the case of 
exogenous propagation of opinion A is studied based on the model with 

Ta = 8 time units , Tg = oo. (34) 

In Figure [7| are presented the simulation results for the parameter values and initial conditions 

n A = 0.3, k b = 0.6, /m(0) = 0.8, p B (0)=0.2, Pn {0) = 0. (35) 



12 



If there were no external forces the trajectory would converge towards the upper left corner cor- 
responding to decision B, as can be seen in the corresponding time-evolution shown in the left 
sub-figure of Figure[71 and in the projection onto the (p a , /Ob)-plane shown in the right sub-figure by 
the thick blue line. For the case that A is propagated through an exogenous force with arbitrary 
18% of the population being impacted each Ta — 8 time units the corresponding time evolution is 
shown in the central sub- figure of Figure [7J and the projection onto the (p a , pb)-plane is shown in 
the right sub-figure by the black line. It can be seen that at each period there is a force pushing 
18% of the population towards decision A. Accordingly, the trajectory does not reach the attractor 
(corresponding to the decision B), but is maintained in an oscillatory regime about some intermedi- 
ate decision distribution. This illustrates the impact of exogenous forces on the dynamic behavior 
in the network and the influence on the overall (mean) decision vector p. 




Time t Time t 



Figure 7: Left: Time evolution of the unperturbed decision dynamics on a dynamic network for the 
parameter values and initial conditions given in (|35[) and a reactive transmission process. Center: 
Timc-cvolution of the associated decision dynamics with exogenous influence on 18% of the popu- 
lation each Ta = 8 time units. Right: Projection onto the (p a , pb)-plane of the trajectories of the 
mean decision probability p (|16p for the unperturbed dynamics (thick blue line), and the one with 
exogenous perturbation (black line). 



3.4 Dynamic network with mass media influence of both opinions 

In order to illustrate the competition of two external forces, consider the case that both decisions 
are propagated through external forces into the network (e.g. by exogenous mass media) but with 
different periods 

Ta = 8 time units, Tb = 16 time units. (36) 

It is considered that both mechanisms affect 10% of the total population to their favor. The 
parameter values and initial conditions were set to 

k a = k b =0-5 7 pi(0) = [0.25,0.15,0.6]', p 2 (0) = [0.15, 0.25, 0.6]' (37) 
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0.2 0.4 0.6 0.8 1 

Pa 

Figure 8: Projection onto the (p a , pb)-plane of two representative trajectories of the mean decision 
probability p (|16[) for the parameter values and initial conditions given in (|37[) and a reactive 
transmission process. The unperturbed trajectories arc represented by the thick blue lines, and the 
ones with exogenous perturbation by the black lines. 

corresponding to two representative cases as shown in Figure |U As illustrated by the thick blue 
lines, in the unperturbed case, the trajectory starting at pi(0) Q37[) converges towards the attractor 
at the upper left corner, and the trajectory starting at /32(C)) (]37|) converges towards the attractor 
at the lower rightt corner. In presence of the exogenous perturbations with periods given in (|36j) 
and affecting both a total number of 10% of the population to their favor, the trajectories are given 
by the black lines. It can be seen that the trajectory starting at /3i(0) does no longer converge 
towards the attractor at the upper left corner, but enters into an oscillatory regime about some 
intermediate mean decision value p, and the trajectory starting at /?2(0) almost converges towards 
the lower right corner but enters into an oscillatory regime close to the corresponding attractor for 
the unperturbed trajectory. 

These results illustrate the differences that may be caused by exogenous propagation mechanisms 
on the decision dynamics in complex networks. In particular it reveals that the final behavior still 
represents nonlinear character implying that the final decision state will strongly depend on the 
initial condition. 

4 Concluding remarks 

The decision dynamics in a scale-free network has been analyzed with respect to inherent nonlin- 
ear convincing mechanisms, impact of network variations over time, and mass media favoritism of 
one single or both opinions. A mathematical model was derived on the basis of some assumptions 
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establishing a quantitative means to analyze the main mechanisms present in decision dynamics. 
Based on a preliminary analysis of the system's fixed points, it was shown that for a static network 
there are between three and four fixed points, as a consequence of the system nonlinearity. Numer- 
ical simulation studies were presented illustrating the predicted behavior for static and dynamics 
scale-free networks, which is in correspondence to the one expected in real networks. The impact 
that mass media have on the decision dynamics was illustrated in different characteristic situations. 
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